Thermodynamics of polymer adsorption to a flexible membrane 



Steffen Karalus/' 2' Wolfhard Janke/'Q and Michael Bachmann^'^^jt] 

^Institut fiir Theoretische Physik, Universitdt Leipzig, Postfach 100 920, D-04009 Leipzig, 
and Centre for Theoretical Sciences (NTZ), Emil-Fuchs-Strafie 1, D-04105 Leipzig, Germany 
^Institut fiir Festkorperforschung (IFF-2) and Institute for Advanced Simulation (IAS-2), 
Forschungszentrum Jiilich, D- 524-25 Jiilich, Germany 
^ Institut fiir Theoretische Physik, Universitdt zu Koln, Zillpicher Strafie 77, D-50937 Koln, Germany 
"^G enter for Simulational Physics, The University of Georgia, Athens, Georgia 30602, USA 

We analyze the structural behavior of a single polymer chain grafted to an attractive, flexible 
surface. Our model is composed of a coarse-grained bead-and-spring polymer and a tethered mem- 
brane. By means of extensive parallel tempering Monte Carlo simulations it is shown that the 
system exhibits a rich phase behavior ranging from highly ordered, compact to extended random 
coil structures and from desorbed to completely adsorbed or even partially embedded conforma- 
tions. These findings are summarized in a pseudophase diagram indicating the predominant class 
of conformations as a function of the external parameters temperature and polymer-membrane in- 
teraction strength. By comparison with adsorption to a stiff membrane surface it is shown that the 
flexibility of the membrane gives rise to qualitatively new behavior such as stretching of adsorbed 
conformations. 



I. INTRODUCTION 

The interaction of macromolecules and cell membranes 
is essential for almost all biological processes. Membrane 
proteins like glycoproteins and transmembrane proteins 
govern the exchange of signals, small molecules, and 
ions between the intracellular and extracellular solvent. 
Membrane-embedded receptors are specific for the bind- 
ing of ligands. The conformational changes caused by the 
binding process can trigger cellular motion, drug delivery, 
or enzymatic catalysis. 

It is interesting to understand how the conformational 
changes that a polymer can experience in the binding 
process to a membrane-like substrate depend on external 
parameters such as temperature and adsorption strength. 
Our study aims at the systematic investigation of such 
transitions for classes of polymer-membrane systems. 
Therefore, it is our goal to construct a conformational 
phase diagram that comprises the generic phase behav- 
ior of these types of organic hybrid systems. 

Initiated by a different motivation, the study of bind- 
ing affinity and specificity of organic and inorganic mat- 
ter, much work has been dedicated to the identification 
of structural transitions polymers experience when ad- 
sorbing to solid substrates [THE]- In these studies, the 
substrate is typically considered as a solid material (e.g., 
a crystal) with virtually no thermal activity (i.e., its sur- 
face structure does not change in the course of thermal 
fluctuations) . 

Although being equally important, much less is known 
about the thermodynamic structural behavior of a poly- 
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mer that interacts with a fluctuating surface such as a 
membrane [13-16 . 

In this paper, we are going to study a simple coarse- 
grained model system consisting of a flexible, elastic 
polymer grafted to a fluctuating substrate by means of 
generalized-ensemble Monte Carlo computer simulations. 
In a comparative analysis, we show how the structural 
phase diagram changes under the influence of thermal 
membrane fluctuations. 

The paper is organized as follows. In Sec. [ll] the model 
system is described in detail. Also, the parallel tempering 
Monte Carlo method is reviewed shortly and the mea- 



sured observables are introduced. Section [III| presents 
and discusses the main results, the pseudophase diagrams 
for the two systems under comparison. Finally, Sec. [TV| 
concludes the paper with a summary of our flndings. 



II. MODEL AND METHODS 

A. Compound polymer— membrane model 

We employ a coarse-grained off-lattice model for a 
single elastic flexible homopolymer |T7H2T| consisting of 
= 13 monomers [22 . All monomers interact pairwise 
via a Lennard- Jones (LJ) potential. 



^LTM = 4epp {a/r)^^-{a/T) 



(1) 



modeling van der Waals forces. Here, r denotes the rel- 
ative distance between the monomers, a = ro/2^/^ the 
zero of the LJ potential, and ro the minimum-potential 
distance. The interaction strength is determined by the 
energy parameter epp. Throughout this study we set 
ro = 1 and epp = 1 as basic length and energy scales. 
Additionally, adjacent monomers are tied together by 
the flnitely extensible nonlinear elastic (FENE) poten- 
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^fene(^) 



{l-[(r-ro)/i?p]'} 



(2) 



modeling covalent bonds, ro is the same minimum- 
potential distance as above. The FENE potential be- 
haves like a harmonic potential with spring constant K 
in the vicinity of ro and diverges for r tq ^ Rp. For 
the simulations we set K = 40 and = 0.3. The con- 
formation of the polymer is completely defined by the set 
of position vectors {Ti}i=i^,,,^N' Eventually, the contri- 
bution of the polymer to the total energy reads 



N-1 N 
i=l j 



N-1 



via 



i+1 



1). 



i+1 



(3) 

The fluctuating substrate is modeled by a tethered 
membrane [25-27 where the individual building seg- 
ments (nodes) are tied together according to a square 
lattice structure with x Ly nodes in total. As tether- 
ing potential we apply, again, a FENE potential. 



FENE 



(r) 



K 



Rl\n[l-[(r-ro)/R^f} . (4) 



The equilibrium distance ro and spring constant K take 
the same values as above, the maximum extension we 
set to = O-l- Additionally, we introduce a hard- 
sphere potential between all pairs of nodes to ensure self- 
avoidance, 



^hs(r) 



if r > 2rhs , 
oo if r < 2rhs , 



(5) 



with hard- sphere radius rhs = 0.15. The configuration 
of the membrane is then described by the set of position 
vectors {r/e, J(/e,^)=(i,i),...,(L,,Lj- The contribution of the 
membrane to the total energy of the system thus reads 

/e=l 1=1 

Ly-1 

+ E E ^F^EdrM-rM+il). (6) 

/c=l 1=1 

In our simulations, we set Lx = Ly = 27. The polymer, 
which is anchored at the membrane center, can then take 
on all possible conformations and its fluctuations are not 
limited by the membrane boundaries. 

The interaction between polymer and membrane is 
modeled by another LJ potential between all pairs of 
monomers and membrane nodes. 



(7) 



The interaction strength between polymer and mem- 
brane is determined by the energy parameter epm which 



serves as an external control parameter, with values be- 
tween 0.05 and 1.50. Furthermore, between the first 
monomer and the central node acts the FENE potential 
^fene(^)' ([H' anchoring potential. The 

interaction term is the third contribution to the energy 
of the system, 

N Ly 



i=l k=l 1=1 

+ '^F%E(kl - ^L,/2,Ly/2\) , 



(8) 



yielding a total energy E = £;p°i + E"^^"^ + E'""^ en- 
tering the partition function Z — ^exp(— where 
/3 = Xjk^T is the inverse thermal energy and the summa- 
tion extends over all possible microstates of the polymer 
and membrane. 

For all the simulations presented in this work, we con- 
strain the membrane to fixed boundary conditions, that 
is, all nodes at the boundary of the membrane are kept 
fixed. Their positions are uniformly arranged on a rect- 
angle of dimensions [L^ — l)ro x [Ly — l)ro in the xy plane, 
such that the membrane can adopt its ground state of 
a regular square lattice. This is a simple choice which 
satisfies computational needs. A different composition 
of the membrane mesh structure will doubtlessly have 
an influence on the location and specific properties of 
low-temperature crystalline pseudophases of the system. 
Nonetheless, we expect that the qualitative structure of 
the pseudophase diagram, as obtained in our study, re- 
mains widely unchanged. The main feature of our mem- 
brane model is its flexibility which allows for thermal 
fluctuations around the planar ground state. To iden- 
tify the effects caused by this flexibility we also study 
the case where the whole membrane is kept fixed in its 
ground state. We will distinguish these two situations 
by speaking of the "flexible membrane" and "stiff mem- 
brane" systems. 



B. Simulation method 

The parallel tempering (PT) Monte Carlo algorithm 
[28l l29j provides a conceptually simple but efficient 
method for the simulation of complex systems over a 
broad temperature range. The basic idea is to perform 
parallel simulations on a number of identical copies (repli- 
cas) of the system at different temperatures using a stan- 
dard sampling scheme (e.g., the Metropolis algorithm 
[30]). After a certain number of sweeps (Monte Carlo 
"time" steps) the conformations are exchanged with an 
acceptance probability 



ipT 



1 for A < 

exp(-A) forA>0, 



(9) 



where A = (/^^ - /3^)[£;(X) - £;(X0]. X and X' are the 
conformations of the replicas at inverse thermal energy 
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Prn and Pn-, respectively. If the temperature values of the 
PT siniulation are chosen thoroughly [31], the conforma- 
tions will perform a random walk in temperature space 
and therefore rapidly decorrelate, helping to escape from 
trapped states. 

The density of states and canonical expectation val- 
ues for temperatures in between the selected values can 
be obtained by the multiple-histogram reweighting tech- 
nique [32]. 



C. Observables 

To identify and distinguish the different structural 
phases of the system we define a number of observ- 
able quantities O. The behavior of the canonical ex- 
pectation values {O) and their temperature derivatives 
d{0)/dT = k^p'^[{OE) - {0){E)] wih provide informa- 
tion about the structural phases and transitions. 

The most basic quantities are the canonical averages 
of the total energy {E) and its individual contributions 
(EP^^), and (^^"*) together with the associ- 

ated heat capacities (temperature derivatives) C, 
C"^^"^, and C^^S where Cp°^ = d {E"^""^) / dT , and so 
on. Evidently, {E) = (£;p°1) + {E"^^"^) + {E'''^) and 

^ ^pol _|_ ^mem _|_ ^int 

The radius of gyration Rg is a measure for the over- 
all compactness of the polymer conformation. It is de- 
fined as the root-mean-square distance of the individual 
monomers from the center of mass of the polymer Tc m = 

N-'j:ti^n, (J?g) = ([iV-^Eti(rn-rc.„..)2]i/2). Al- 
though the membrane surface is not planar, it may be in- 
structive to separate the radius of gyration into the com- 
ponents of the gyration tensor parallel and perpendicular 
to the membrane equilibrium state, i.e., to the x?/-plane, 
= {{N-^Y.n=i{i^n-x,.^.f + {y^-y,.^.f]r/^) 

and (i?g,i) = {[N-^En=ii^n - ^c.m.)']'/'}. The ratio 
of these two values gives the sphericity aspect ratio with 
respect to the x?/-plane = V^{Rg,i.)/{Rg,\\)' This def- 
inition is chosen such that = 1 indicates spherically 
symmetric structures, while an oblate (prolate) spheroid 
with polar axis z will have < 1 (^r > !)• 

To determine whether the polymer is, on average, close 
to the membrane surface (adsorbed) or freely exploring 
the third dimension (desorbed), we measure the distance 
of the center of mass of the polymer from the membrane 
equilibrium state in the xy plane (i.e., its z component. 

The number of contacts between monomers and 
membrane nodes may give excellent information about 
the state of adsorption. As a reasonable, but 
still to some extent arbitrary, measure for a "con- 
tact" we decided to count every monomer-node 
pair that contributes a V^^{r) energy less than 
a threshold value E^^ = — 0.5epm as a con- 
tact yielding a mean number of polymer-membrane 

contacts (np^) = (^-^ E^=i E^ii Efji ©(^r - 



^Lj^d^* ~^k,i\)))^ where 9(x) denotes the Heaviside step 
function. Similarly, we define the number of intrin- 
sic contacts as another measure for the compactness 
of the polymer, (npp) = (iV"^ E^T' Ef=.+i ©[^c^^ - 
^LT(|rz-r,l)]), with ^PP^ -0.56pp. 



III. RESULTS AND DISCUSSION 
A. Stiff membrane system 

The stiff membrane system is studied as a reference 
system here. On the one hand, the knowledge about 
this system will enable us to point out the differences in 
the behavior of the stiff and the fiexible membrane cases 
and to identify the new effects emerging from the surface 
flexibility. On the other hand, the stiff membrane case 
allows a comparison of our model with studies on poly- 
mer adsorption to solid (flat) substrates |TT1[T2]. For the 
stiff membrane system we keep all membrane nodes, in 
agreement with the boundary conditions, flxed on a reg- 
ular square lattice with lattice spacing tq in the xy plane. 
This is obviously the ground state of the membrane as 
all FENE springs are at equilibrium, "I^^neI^o) = 0, 
^mem _ deviation from this state would in- 

troduce a positive contribution to the membrane energy. 

The main results we show here were obtained by par- 
allel tempering simulations with 24 replicas in the tem- 
perature range from 0.021 to 1.500 and 8 x 10^ sweeps 
on each replica. Exchanges of conformations between the 
replica were attempted every 20 sweeps. For the uncriti- 
cal region above T = 1.5, additional PT simulations with 
10^ sweeps were carried out. All the simulations were 
performed at 30 different values of epm in the interval 
ep^ = 0.05,..., 1.50. 



1. Pseudophase diagram 

The main information about the structural behavior is 
summarized in the pseudophase diagram shown in Fig. [l] 
It displays the structural pseudophases and pseudo-phase 
transitions in the Cpm-^-plane. Temperature increases 
from bottom to top, interaction strength increases from 
left to right. 

The structural phases in which certain conformations 
predominate are labeled by a letter code adopted from 
Ref. [11 . Since we are working with a flnite system, we 
cannot identify precise transition lines but rather tran- 
sition regions between the phases. These are displayed 
as the shaded regions in the phase diagram where the 
dark shades indicate well-founded transitions and the 
bright shades stand for less established transitions. In the 
following, we describe all the identifled structural pseu- 
dophases, for pictures of typical conformations in each 
phase, see Fig. |2j 
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FIG. 1. Pseudophase diagram for the stiff membrane sys- 
tem. The main structural phases divide into predominantly 
adsorbed (A), and desorbed (D) on the one hand and ex- 
panded (E), globular (G), and compact (C) structures on the 
other hand. 




(a) DC 



(b)G 



(c) DE 




(d) ACl 



(e) AC2 



(f) AE2 




(g) AEl 

FIG. 2. (Color online) Representative conformations in the 
individual phases of the stiff membrane system (see text for 
a detailed description). 



Desorbed Compact (DC): The most compact confor- 
mations the polymer can adopt. For the exam- 
ple considered here {N = 13), these are the highly 
ordered icosahedral states we know from the free 
polymer [19|. Here, it is simply attached to the 
membrane surface [Fig. |2|a)]. 



Globular (G): Compact but disordered conformations. 
Similar to DC these are globular states of the free 
polymer attached to the membrane. Because of 
the short chain length, we cannot distinguish here 
between adsorbed and desorbed globular conforma- 
tions [Fig. ^ 



Desorbed Expanded (DE): The most disordered con- 
formations. These are random coil structures of 
the free polymer anchored to the membrane and 
restricted to the upper half-space [Fig. [2|c)]. 

Adsorbed Compact, Single Layer (ACl): 

Completely adsorbed, compact conformations. 
These film-like structures are constrained to lie flat 
on the surface of the membrane and to its square 
lattice structure. Additionally they are mainly 
disk-like compact [Fig. [2|d)]. 

Adsorbed Compact, Double Layer (AC2): 

Partially adsorbed, compact conformations. 
These consist of a bottom layer adsorbed to the 
membrane and adapted to its lattice structure and 
a second layer on top of it building a semi-spherical 
droplet [Fig.[2];e)]. 

Adsorbed Expanded, Double Layer (AE2): 

Partially adsorbed, extended conformations. 
These are random coil structures mainly adsorbed 
to the membrane surface but also extending into 
the third dimension [Fig. [2| 



Adsorbed Expanded, Single Layer (AEl): 

Completely adsorbed, extended conformations. 
These are mainly two-dimensional random coil 
structures adsorbed to the membrane surface 
[Fig.gg)]. 



2. Structural phases and transitions 

In the following, we substantiate the proposed pseu- 
dophases and discuss the transitions by looking in detail 
at the observables measured in the course of the simula- 
tions. 

The simulation methods we apply here are not spe- 
cialized in finding ground states. For this task simpler 
and more efficient algorithms exist. Nevertheless, our 
simulation results reach down to sufficiently low temper- 
atures where a convergence sets in. We therefore think 
that the low-energy conformations we found in the low- 
temperature phases DC, AC2, and ACl are good approx- 
imations of the real ground states. Of course, further de- 
tails like the actual bond distribution, translations on the 
membrane surface, and distortions of the idealized sym- 
metries due to LJ interactions over longer ranges have to 
be considered for more precise ground-state predictions. 

The different compact conformations are characterized 
in a simple way by the contact numbers (see Sec. II C). 
The icosahedron is the most compact structure of the 



5 



ft 
ft 



ft 



0.1 



0.4 ^0.7' 



1.0 ^ 1.3 

0.2 ^0.5 0.8 ^1.1 ^1.4 
0.3 ^0.6 0.9 ^1.2 ^1.5 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.^ 

T 




ft 



70 
60 
50 
40 
30 
20 
10 



' 1 ' 


1 ■ 1 ■ 1 ■ 

free polymer 


• 










\//\\ 0.10 




/ 


a/ \ \ 0-20 


□ 




n \\ 0-30 


■ 






o 


1 


I.I.I. 





^Ife 0.4 \- 



0.2 - 




FIG. 3. (Color online) Mean number of intrinsic (top) and 
polymer-membrane (bottom) contacts as a function of tem- 
perature for various polymer-membrane interaction strengths 
in the stiff membrane system. The legend applies to both 
plots. 



FIG. 4. (Color online) Polymer heat capacity (top) and tem- 
perature derivative of the radius of gyration (bottom) for low 
polymer-membrane interaction strengths in the stiff mem- 
brane system. For better comparison, the values for the free 
polymer are also plotted (red circles). 



13mer in DC, resulting in the highest number of intrinsic 
contacts. Counting them correctly we end up with n^^ = 
42/13 ~ 3.23. Due to its compactness, this structure can 
build up few polymer- membrane contacts, ~ 1. 

In the two-dimensional compact disk-like structures 
of ACl every monomer is in contact with four mem- 
brane nodes, = 4, and we count n^^^ = 18/13 ^ 
1.38 intrinsic contacts per monomer. In the partially 
adsorbed semi- spherical droplets of AC2, nine of the 
monomers are in contact with four membrane nodes each, 
^pm^ — 36/13 ^ 2.78, for the intrinsic contacts we get 
Tipp^ = 32/13 ~ 2.46. Equipped with this informa- 
tion, we can easily identify the different low-temperature 
phases. Figure |3] shows the mean numbers of intrinsic 
and polymer-membrane contacts for different polymer- 
membrane interaction strengths over the whole interac- 
tion strength spectrum investigated here. At low temper- 
atures, we see convergence to exactly the values predicted 
above. 

For epm < 0.3, the mean numbers of intrinsic contacts 
approach the value 3.23 uniquely characterizing DC. The 
mean number of polymer-membrane contacts does not 



exhibit such a clear convergence in this case, since re- 
orientation on top of the membrane is still possible in 
this phase. For 0.4 < < 0.8, the mean contact num- 
bers approach 2.46 for the intrinsic contacts and 2.78 for 
the polymer-membrane contacts. Together, these values 
identify compact double-layer conformations as approxi- 
mate ground states in this phase (AC2). For epm > 0.9, 
the mean numbers of 1.38 intrinsic contacts and four 
polymer-membrane contacts are approached for low tem- 
peratures. We therefore find monolayer structures as ap- 
proximate ground states in this region (ACl). 

For low polymer-membrane interaction strength, the 
polymer should behave similarly to the free polymer 
p!8| ITQ] . The freezing transition between the compact 
icosahedral and the globular phase is identified by peaks 
in the polymer heat capacity and temperature derivative 
of the radius of gyration around T = 0.35. The latter 
quantity shows another peak at around T = 0.9 indicat- 
ing the 6 transition between the globular and expanded 
phase. In Fig. [4} we plot these two quantities for the stiff 
membrane system. The curves for epm = 0.05 and 0.10 
lie almost on top of each other and agree almost perfectly 
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FIG. 5. (Color online) Contribution of polymer- membrane 
interaction to heat capacity (top) and sphericity aspect ratio 
(bottom) at the AC2 - G transition in the stiff membrane 
system. 
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FIG. 6. (Color online) Temperature derivatives of the parallel 
(upper curves) and perpendicular (lower curves) component 
of the gyration tensor at the AC2 - AE2 transition in the stiff 
membrane system. 
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with those for the free polymer (red circles). We conclude 
that the same kinds of transitions take place here, from 
DC over G to DE. For = 0.20 and 0.30 the freezing 
peaks shift to lower temperatures in both quantities, dis- 
appearing almost completely at epm = 0.35 which marks 
the point where the ground state changes from DC to AC2 
(see Fig. [3|. The 9 peak maintains its position T ^ 0.9 
in the whole range of 0.05 < epm < 0.35. 

In the range of 0.4 < epm ^ 0.6, the mean polymer- 
membrane contacts decrease considerably between T ^ 
0.15 and 0.35 without a substantial change in the in- 
trinsic contacts (Fig. |3| suggesting a (partial) desorption 
transition into the globular phase G. This assumption is 
supported by a peak in the interaction part of the heat 
capacity C^"*(T), indicating that the interaction energy 
becomes higher in this region, and an increase of the 
sphericity aspect ratio in that range (Fig. [5|. We find 
that the transition temperature becomes higher for in- 
creasing interaction strength, from T ^ 0.18 at epm = 0.4 
toT ^ 0.30 at 

^pm — 0.6. For larger values of the interac- 
tion strength, we see the sphericity aspect ratio decreas- 
ing which suggests a transition into structures extending 
parallel to the membrane. 



FIG. 7. (Color online) Temperature derivative of the parallel 
component of the gyration tensor at the G - AE2 transition 
in the stiff membrane system. 



This can be affirmed by looking at the two compo- 
nents of the gyration tensor separately. Figure [6] displays 
the temperature derivatives of the two components for 
0.7 < 

^pm ^0.8. We see a clear peak in the parallel com- 
ponent (upper three curves) and no signal or a valley in 
the perpendicular component (lower three curves) around 
0.30 ^ T < 0.35 meaning that the polymer does extend 
in the xy direction while it does not, or even contracts, in 
the z direction. This strongly supports the assumption 
that we have a transition from the droplet phase AC2 to 
the adsorbed three-dimensional random coil phase AE2 
here. 

The transition from the globular phase G to the ad- 
sorbed expanded, double layer, phase AE2 is difficult to 
identify since globular structures are due to the anchor- 
ing approximately as close to the surface as the adsorbed 
but still three-dimensional random coils of AE2. Never- 
theless, we do observe an increase in the parallel com- 
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FIG. 8. (Color online) Temperature derivative of the z compo- 
nent of the center of mass (top) and perpendicular component 
of the gyration tensor (bottom) at the AE2 - DE transition in 
the stiff membrane system. 



ponent of the gyration tensor indicated by a peak in its 
temperature derivative for epm ^ 0.6 (Fig. [t]). For in- 
creasing epm, the peak position slightly moves to lower 
temperatures and disappears from epm ^ 0.7 on. We 
will see below that the perpendicular component experi- 
ences a similar increase, which for 0.3 < epm ^0.6 takes 
place at significantly higher temperatures. Therefore, we 
conclude that the G - DE transition discussed above sepa- 
rates for higher interaction strengths into two transitions, 
an expansion G - AE2 and a desorption AE2 - DE. 

The suggested desorption transition from adsorbed 
three-dimensional random coils AE2 to desorbed random 
coils DE should be measurable by a simultaneous increase 
in the perpendicular component of the gyration tensor 
and the height of the polymer's center of mass above the 
membrane surface. These two effects can indeed be ob- 
served as peaks in the respective temperature derivatives 
(Fig. [8|. As in many transitions of finite systems, the 
peak positions slightly differ from each other causing the 
rather broad transition region in the phase diagram. 

We conclude this section by discussing the last tran- 
sition between the compact film phase ACl and the ad- 
sorbed expanded phases AE2 and AEl. Figure [9] shows 




FIG. 9. (Color online) Interaction heat capacity at the ACl 
AE2 transition in the stiff membrane system. 



the interaction heat capacity for high polymer-membrane 
interaction strengths 1.0 < epm < 1-5. We see a clear 
peak indicating a (partial) desorption transition shifting 
to higher temperatures for stronger polymer-membrane 
interaction. This is supported by peaks in the tempera- 
ture derivatives of the components of the gyration ten- 
sor (Fig. 10). Taking a closer look at these figures, we 
see that with increasing interaction strength the paral- 
lel expansion shifts more slowly to higher temperatures 
than the expansion perpendicular to the membrane. Al- 
though this disagreement might be another finite-size ef- 
fect, we think that it is caused by the onset of the com- 
pletely adsorbed random coil phase AEl. The existence 
of this phase is also supported by the fact that the paral- 
lel expansion becomes stronger while the perpendicular 
expansion becomes weaker for higher polymer interaction 
strength, as suggested by the peak heights. 



3. Comparison with established results 

The phase diagram we have constructed here (Fig. [l]) 
qualitatively compares with that by Moddel et al. [llj. 
Fig. 2, which allowed us to adopt their notation. The 
main qualitative differences between the models intro- 
duced there and here are the anchoring of the polymer 
and the discrete surface structure. For this reason, we do 
not expect quantitative agreement in the locations of the 
transitions but qualitatively, both systems should behave 
similarly. 

As mentioned above we cannot distinguish between ad- 
sorbed and desorbed globular structures in our case, so 
we just have a single globular phase G in contrast to the 
adsorbed globular AG and desorbed globular DG phases 
in Ref. [llj. Also we do not find an adsorbed compact 
phase, like AC2a in Ref. [11], that can be clearly distin- 
guished from DC and AC2. Obviously these two differ- 
ences stem from the anchoring of the polymer. For longer 
chains the distinctions possibly can be recovered also in 
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FIG. 10. (Color online) Temperature derivative of the parallel 
(top) and perpendicular (bottom) components of the gyration 
tensor at the ACl - AE2 / AEl transition in the stiff membrane 
system. 



the grafted case. 

Another difference can be identified in the transition 
line between compact and globular or expanded struc- 
tures. In Ref. [11], these transitions occured at more or 
less constant temperature, independently of the polymer- 
surface interaction strength. This seems to be not the 
case in our model. 



B. Flexible membrane system 

We now extend the focus of this study and discuss the 
structural behavior of the coupled system consisting of an 
elastic polymer anchored to a flexible membrane. This 
section describes our findings about the flexible mem- 
brane system in detail where we concentrate on the dif- 
ferences and similarities with the stiff membrane case. 

All results presented here were obtained by parallel 
tempering simulations with 48 replicas in the tempera- 
ture range from 0.055 to 1.500 and 10 sweeps on each 
replica. Exchanges of conformations between the replica 
were attempted every 20 sweeps. The simulations were 




FIG. 11. Pseudophase diagram for the flexible membrane sys- 
tem. The main structural phases divide into predominantly 
adsorbed (A), desorbed (D), and embedded (M) on the one 
hand and expanded or elongated (E), globular (G), and com- 
pact (C) structures on the other hand. 



As the performance of the parallel tempering conforma- 
tion exchanges turned out to be poor in the temperature 
range below T ^ 0.1 leading to large statistical errors, we 
restrict our discussion to temperatures above this value. 



1. Pseudophase diagram 



carried out at 30 different values of er. 



0.05,..., 1.50. 



The pseudophase diagram shown in Fig. 11 summa- 
rizes the main information about the structural behavior 
of the flexible membrane system. It displays the struc- 
tural pseudophases and pseudo-phase transitions in the 
^pm-T plane. As before, temperature increases from bot- 
tom to top and polymer-membrane interaction strength 
from left to right. 

Again, the structural phases are labeled by a letter 
code similar to the one used before and transition re- 
gions are shaded in dark gray for well-founded and light 
gray for less certain transitions. In the following, we 
describe the identified pseudophases in the flexible mem- 
brane system. Images of typical conformations in each 
phase are shown in Fig. [12] The phases in which the 
polymer-membrane interaction is not of particular im- 
portance (DC, G, and DE [Fig. [T2|a)-(c)]) are, of course, 
very similar to the case of the stiff membrane. To avoid 
redundancy we omit a description of these here. The 
main new structures are the embedded conformations 
(MC, ME) which reflect the back-reaction between poly- 
mer and membrane fluctuations. Let us now first sum- 
marize the differences in the phase behavior, compared 
to the stiff membrane case. 

Adsorbed Compact, Single Layer (ACl): 

Completely adsorbed, compact conformations. 
These film-like structures lie flat on the membrane 
surface. Since the membrane is deformable in this 
case, the conformations are, strictly speaking, not 
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(d) ACl (e) MC (f) ME 




(g) AEl (h) AE2 

FIG. 12. (Color online) Representative conformations in the 
individual phases of the flexible membrane system (see text 
for a detailed description). 



two-dimensional as before and do not necessarily 
adopt any lattice structure [Fig. [T2|d)]. 

Embedded Compact (MC): Compact conformations, 
deforming the membrane. These highly com- 
pact structures are partially wrapped in by the 
membrane which, on the other hand, forms an 
immersion to incorporate the collapsed polymer 
[Fig. [He)]. 

Embedded Elongated (ME): Stretched conforma- 
tions, incorporated by the membrane. For these 
structures the membrane forms a channel into 
which the almost fully stretched polymer is em- 
bedded. The polymer is not randomly expanded 
here, but specifically elongated [Fig. [l2|f)]. 

Adsorbed Expanded, Single Layer (AEl): 

Completely adsorbed, extended conformations. 
These are randomly expanded structures that live 
on the membrane surface. In contrast to the stiff 
membrane case, the structures are not strictly 
two-dimensional as the membrane is fluctuating in 
all directions [Fig. [l2|g)]. 

Adsorbed Expanded, Double Layer (AE2): 

Partially adsorbed, extended conformations. 
These are partly adsorbed random coil structures 
significantly extending into the third dimension. 
Although the membrane is important here, there 
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FIG. 13. (Color online) Polymer heat capacity (top) and tem- 
perature derivative of the radius of gyration (bottom) for low 
polymer-membrane interaction strengths in the flexible mem- 
brane system. For better comparison, the values for the stiff 
membrane system at epm = 0.05 (Fig. |4]) are plotted again 
(red circles). The legend applies to both plots. 



is little structural difference to the stiff membrane 
case [Fig.[l2];h)]. 



2. Structural phases and transitions 

We will shortly summarize similarities with the stiff 
membrane system here and then look at the differences 
and emergence of new phenomena in more detail. 

The phases DC, G, and DE are hardly affected by 
the presence of the membrane. Therefore, we expect 
the fluctuations of the membrane to have little effect 
on these phases and the transitions between them. For 
the adsorbed random coil phase AE2, the membrane is, 
of course, essential as the adsorbing surface. Neverthe- 
less, the phase is already highly disordered such that 
membrane fluctuations should not generate new effects 
here. This is exactly what we observe in our simula- 
tion results and is reflected in the phase diagram where 
we just have marginal differences to the stiff membrane 
system. We do not repeat all the plots of the observ- 
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FIG. 14. (Color online) Mean number of intrinsic (top) and 
polymer-membrane (bottom) contacts as a function of tem- 
perature for various polymer-membrane interaction strengths 
in the flexible membrane system. The legend applies to both 
plots. 



ables identifying the phases just mentioned and their 
transitions. To provide only one example for these obvi- 
ous similarities, Fig. [l3] shows the polymer heat capacity 
and the temperature derivative of the radius of gyration 
for low interaction strengths, 0.05 < Cpm < 0.45. The 
curves clearly resemble those in Fig. |4] for the stiff mem- 
brane system indicating, again, the O collapse around 
T ^ 0.9 and the freezing transition at T ^ 0.3 of the 
polymer. Note that we already find a small difference 
to the stiff membrane system here. The freezing peak 
is still pronounced for Cpm = 0.35 and present up to 
Cpm = 0.40 where this transition already disappeared 
in the stiff membrane system. This possibly indicates 
a slightly higher stability of the icosahedral phase in the 
flexible membrane system. A reason for this could be the 
membrane's possibility to adapt to the polymer structure 
and build up more polymer-membrane contacts while the 
polymer maintains its icosahedral shape. 

In the stiff membrane system, the mean numbers of in- 
trinsic and polymer-membrane contacts showed a strik- 
ing convergence to our predicted values for the highly 
ordered low-temperature states. For the case of a flex- 
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FIG. 15. (Color online) Sphericity aspect ratio at the MC - G 
and MC - ACl transitions in the flexible membrane system. 



ible membrane, the curves do not exhibit such a clear 
convergence (Fig. [l4|, but still, we can extract some in- 
formation about the low-temperature behavior of the sys- 
tem. For Cpm < 0.3, again we observe convergence to a 
mean of 3.23 intrinsic contacts uniquely characterizing 
the icosahedral phase DC. 

For epm > 1.0, the mean number of polymer- membrane 
contacts increases significantly above 4 while the mean 
number of intrinsic contacts drops down to values even 
below those of the disordered random coil phases. The 
latter number suggests that the probability of specifically 
stretched structures is enhanced in this region, whereas 
in random coils expanded structures only dominate be- 
cause of higher entropy. The former number can only be 
reached when the membrane locally wraps around the in- 
dividual monomers allowing for more than four polymer- 
membrane contacts. In the most perfect realization of a 
structure in the embedded elongated phase ME, the poly- 
mer would be linearly stretched, yielding 12/13 ^ 0.92 
intrinsic contacts, and incorporated in a channel with 
a triangular cross section such that every monomer is 
in optimal distance to two square patches of the mem- 
brane, resulting in a number of six polymer-membrane 
contacts per monomer. Obviously, this idealized situa- 
tion is not reached here but from the tendencies shown 
in Fig. [M] we propose embedded elongated conformations 
(ME) as predominant states for low temperature and high 
polymer-membrane interaction strength. In the range of 
0.7 < Cpm ^ 0.8 the two contact numbers seem to con- 
verge to certain intermediate values. Although we are 
not sure which would be the perfect symmetry in this 
case, we expect to find the embedded compact phase MC 
as the low-temperature phase here as suggested by the 
snapshot in Fig, p^e). 

The transition out of the MC phase with increasing 
temperature is analyzed in Figs. p!5}|T7l In the range 
of 0.65 < 6pm ^ 0.75, the perpendicular component of 
the gyration tensor significantly increases in the region 
0.14 < T < 0.18 as indicated by the peaks in the cor- 
responding temperature derivative (Fig. 16). Addition- 
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FIG. 16. (Color online) Temperature derivative of the per- 
pendicular component of the gyration tensor at the MC - G 
transition in the flexible membrane system. 
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FIG. 18. (Color online) Temperature derivatives of the radius 
of gyration (top) and the mean number of polymer-membrane 
contacts (bottom) at the ME - ACl transition in the flexible 
membrane system. 



FIG. 17. (Color online) Temperature derivative of the parallel 
component of the gyration tensor at the MC - ACl transition 
in the flexible membrane system. 



ally, we observe a slight increase in the sphericity aspect 
ratio as depicted in Fig. [15] meaning that conformations 
change from oblate to more spherical shapes. This marks 
the transition from embedded compact (MC) to globular 
(G) conformations which qualitatively compares to the 
droplet - globule transition (AC2 - G) in the stiff mem- 
brane case. 

For higher polymer-membrane interaction strengths, 
0.80 < epm ^ 0.95, we observe a peak in the temperature 
derivative of the parallel component of the gyration ten- 
sor in the region 0.15 < T < 0.19 (Fig. [l7|. As the per- 
pendicular component does not increase simultaneously, 
this indicates a further "flattening" of the polymer struc- 
tures. This can also be concluded from the drop in the 
sphericity aspect ratio (Fig. 15). Together with the ob- 
served drop in the number of intrinsic contacts (Fig. 14), 
we identify the transition from the embedded compact 
(MC) to the adsorbed compact (ACl) phase here. 

One of the strongest transition signals we find is shown 



in Fig. 18 The temperature derivatives of the radius 
of gyration and the number of polymer-membrane con- 
tacts exhibit a deep valley indicating a sharp transition 
line. The position of the valley increases from T ^ 0.14 
at epm = 1.0 to T ^ 0.3 at epm = 1.5 almost linearly. 
The decrease in the radius of gyration suggests that we 
have a transition from expanded (or elongated) struc- 
tures at lower temperature to more compact structures 
at higher temperature. This is exactly the opposite of 
what we usually observe in transitions like the O collapse 
and would not be possible without the membrane serv- 
ing as an agent or a medium enhancing the probability 
of elongated conformations in the low-temperature phase. 
Looking also at the mean number of polymer-membrane 
contacts per monomer (Fig.[l4|, we see that the "inverse 
collapse" is accompanied by a drop from more than four 
to slightly less than four contacts per monomer, which 
was the maximum value for a stiff membrane. This sug- 
gests that the membrane adapts its shape to the poly- 
mer to form up to six contacts with each monomer in 
the low-temperature phase and gives up this behavior at 
higher temperatures after passing the transition. Taken 
together, we observe a transition form elongated con- 
formations of the polymer which are partially incorpo- 
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ration tensor at the ACl - AE2 / AEl transition in the flexible 
membrane system. 



rated into the membrane (ME) to disk-like compact film 
structures adsorbed to the fluctuating membrane surface 
(ACl). 

For higher temperatures we observe another, but much 
weaker, transition signal in the temperature derivative of 
the radius of gyration. Taking a closer look at the tem- 
perature derivatives of the individual components of the 
gyration tensor (Fig. 19), we find a similar behavior to 
what we have seen already in the stiff membrane system. 
The perpendicular component increases at a transition 
line ranging from T ^ 0.35 at = 1.0 to T ^ 0.85 



at Cpm = 1.5 indicating a (partial) desorption of the 
polymer from the membrane surface. Additionally, the 
parallel component increases along a transition line from 
T ^ 0.4 at = 1.2 to T ^ 0.55 at epm = 1-5, indi- 
cating an expansion of the polymer parallel to the mem- 
brane surface. As in the stiff membrane case, we con- 
clude that what we observe here is the transition from 
disk-like compact film structures adsorbed to the mem- 
brane surface (ACl) to expanded random coil structures, 
partially adsorbed to the membrane surface (AE2). For 
high polymer- membrane interaction strengths, we think 
that additionally (almost) completely adsorbed random 
coil conformations (AEl) appear in between. 



In the region of the phase diagram where the phases 
G, MC, ACl, and AE2 meet, epm - 0.7, ...,1.0 and 
T ^ 0.2, . . . , 0.6, many transitions seem to superimpose. 
Therefore, it cannot be clearly decided what structures 
are predominantly present in that region. As this is a 
finite-size effect, we think that investigations on larger 
systems could reveal more about the actual phase struc- 
ture here. 



IV. SUMMARY 

In this paper we have studied the structural behavior 
of a coarse-grained model system consisting of a single 
elastic, flexible polymer in contact with a flexible mem- 
brane. The main goal was to identify the conformational 
thermodynamic pseudophases of the system in depen- 
dence of temperature and polymer-membrane interaction 
strength as external parameters. To clearly point out the 
effects emerging from the flexibility and adaptivity of the 
membrane, we separately investigated the two cases of 
a stiff membrane forced into a flat state and a flexible 
membrane, both with the polymer grafted. 

The pseudophase diagrams were constructed from 
the measured observables and their fluctuations for an 
intermediate-sized system of polymer length 13 and 
membrane size 27 x 27. The observed conformational 
pseudophases include the following: 

• Desorbed conformations for low polymer- 
membrane interaction strengths. These phases are 
well-known from studies of free polymers. With 
decreasing temperature the polymer undergoes 
two transitions, the O transition from extended 
random coil structures (DE) to more compact but 
disordered globular conformations (G) and the 
freezing transition to a highly compact icosahedral 
ground state (DC). 

• Adsorbed disordered conformations in both sys- 
tems. In these conformations the polymer is ad- 
sorbed to the membrane surface but still highly 
disordered such that the presence or absence of 
membrane fluctuations does not induce qualitative 
differences. We distinguish between (almost) com- 
pletely adsorbed, extended (AEl) and partially ad- 
sorbed, extended (AE2) random coil structures. 

• Adsorbed well-ordered conformations in the stiff 
membrane system. For low temperatures, we ob- 
served highly ordered conformations adsorbed to 
the membrane surface. We distinguish between 
disk-shaped compact film structures (ACl) and 
semi-spherical droplets (AC2). 

• Embedded conformations in the flexible membrane 
system. At low temperatures, we found the mem- 
brane adapting its structure such that it partially 
incorporates the polymer. Compact oblate-shaped 



13 



structures embedded into the membrane (MC) were 
observed at intermediate and linearly stretched em- 
bedded structures (ME) were identified at high in- 
teraction strength. 

Because of its smallness, the system's structural be- 
havior is influenced by flnite-size effects. Nevertheless, 
we expect that the phase behavior of the system remains 
qualitatively intact also for larger systems. A verifica- 
tion by a systematic finite-size analysis is currently un- 
achievable for most of the phases, in particular in the 
low-temperature regime. Thus, this remains future work 



in its own right. However, also the finite-size effects are 
of substantial interest since classes of polymers, such as 
proteins, are naturally finite. 
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